Gravitational Radiation from Plunging Orbits 
- Perturbative Study - 



Yasushi Minc0 and Jeandrew Brin]<0 
California Institute of Technology, MC 130-33, Pasadena, CA 91125 

Numerical relativity has recently yielded a plethora of results about kicks from spinning mergers 
which has, in turn, vastly increased our knowledge about the spin interactions of black hole systems. 
In this work we use black hole perturbation theory to calculate accurately the gravitational waves 
emanating from the end of the plunging stage of an extreme mass ratio merger in order to further 
understand this phenomenon. This study focuses primarily on spin induced effects with emphasis 
on the maximally spinning limit and the identification of possible causes of generic behavior. 

We find that gravitational waves emitted during the plunging phase exhibit damped oscillatory 
behavior, corresponding to a coherent excitation of quasi-normal modes by the test particle. This 
feature is universal in the sense that the frequencies and damping time do not depend on the orbital 
parameters of the plunging particle. Furthermore, the observed frequencies are distinct from those 
associated with the usual free quasi-normal ringing. Our calculation suggests that a maximum in 
radiated energy and momentum occurs at spin parameters equal to a/M = 0.86 and a/A4 = 0.81, 
respectively for the plunge stage of a polar orbit. The dependence of linear momentum emission on 
the angle at which a polar orbit impacts the horizon is quantified. One of the advantages of the 
perturbation approach adopted here is that insight into the actual mechanism of radiation emission 
and its relationship to black hole ringing is obtained by carefully identifying the dominant terms in 
the expansions used. 



I. INTRODUCTION 

As the field of numerical relativity matures, the gravi- 
tational wave community gains insight into the most dy- 
namic regions of spacetime pj] . Possibly the most spec- 
tacular strong field event in gravity, namely that of a bi- 
nary black hole merger, has recently yielded a series of in- 
teresting results related to the spins of the holes involved 
in the collision. A particular highlight is the discovery 
of the so-called super-kick configuration, where the spin 
interaction drasticallychanges the trajectory of the resul- 
tant black hole 0, 0, 13l ■ The body of knowledge relating 
to spinning compact objects has greatly increased as a re- 
sult, inviting analytic explanations. Important analytic 
strides in this regard have been made [3 H, 0, H, H, [l3| . 

The evolution of a black hole binary proceeds through 
three stages: the inspiral on a quasi stable orbit, the 
plunge and merger and the final ringdown stage. It is in 
the plunge and merger where the strong field fully nonlin- 
ear attraction of general relativity completely dominates 
the dynamics. 

Of these three stages, the first and third have been 
carefully understood analytically. The first stage, during 
which the slow inspiral over a long timescale results in the 
emission of huge amounts of gravitational radiation that 
dominate the observed signal, has been extensively mod- 
eled [m, The final stage, where the highly distorted 
black hole approaches a Kerr black hole via quasi-normal 
ringing, is also carefully understood in terms of vacuum 
black hole perturbation theory [l^, [ij, [l^, [l^ . 
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The dynamics of the intermediary phase of plunge and 
merger, where the binary transitions from a system de- 
scribed by several parameters related to the two black 
holes to a resultant system described by far fewer pa- 
rameters (namely the spin, mass and position of the re- 
sultant single black hole), has been less throughly ex- 
plored. This transition should by its very nature dis- 
play universal behavior. The interest in this phase of 
the binary evolution is threefold. Firstly, the transition 
from inspiral to ringdown marks a turning point in the 
amplitude of the gravitational waveform and is essen- 
tial for the production of accurate templates for obser- 
vation. Secondly, it is thought that the plunging phase 
strongly influences final kick velocity attained by the final 
black hole. Thirdly, this phase gives us the opportunity 
to probe highly nonlinear/non-Newtonian events such as 
horizon formation and the merger process itself, to which 
currently only numerical relativity gives access. 

The first issue can be addressed by attempting to 
extrapolate waveforms produced for the first and third 
stages to bridge the chasm of understanding present in 
the second stage. A method which extends the post- 
Newtonian (PN) approximation into this regime while 
leaning on numerical relativity results was proposed in 
@i 01) however it is mentioned there that the uncer- 
tainty due to the plunging stage remains considerable. 
The amount of linear momentum emitted during the 
plunge stage is expected to be significant in compari- 
son to that emitted during inspiral, making the contri- 
bution from the final plunge stage to the final velocity 
of the kicks of the black hole important. Furthermore, 
this statement also implies that the position in the or- 
bit at which a binary merges strongly influences the end 
result. While the kick velocities resulting from spin in- 
teractions can be modeled effectively by fitting formulas 
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which treat the strong field interactions as an efi^ective 
black box 0, d, 0, 0| , a careful understanding of the in- 
teractions that produce them and the precise connection 
to gravitational radiation still remains somewhat illusive. 
A constructive formalism for treating and understanding 
strongly nonlinear effects in the highly dynamic regions 
of spacetime during the plunge and merger phases is in 
its early stages of development. Subtle details, exhibited 
in the waveform and dynamics during the merger phase, 
such as the anti-kick observed by [l3| give us further clues 
as to dynamics that require explanation. As numerical 
relativity now takes the lead in ushering in a new stage in 
our analytic exploration of the Einstein field equations, 
we can begin to build our intuition about how interac- 
tions between compact object take place and develop the 
tools for describing them. 

In this paper, we study gravitational waves emitted 
from the end of the plunging phase of a small black hole 
falling into a spinning black hole. The purpose of the 
calculation is to obtain a physical understanding of the 
origin of the radiation emitted and to identify possible 
causes of universal behavior. We explore the extreme 
mass ratio case via black hole perturbation theory. This 
approach is used for two reasons. Firstly, it makes the 
problem analytically tractable at minimal computational 
cost. Secondly, it gives insight into the physical mecha- 
nism by which the radiation is produced by the plunging 
phase of the black hole merger and directly relates that 
picture to the resulting gauge invariant radiation quan- 
tities computed at infinity. 

The proposed model allows the effect of spin to be ex- 
plored by considering the background spacetime to be a 
Kerr black hole with mass M and spin a. The second 
smaller black hole is modeled as a point particle mov- 
ing on a geodesic orbit described in SecHIl Past studies 
in which black hole perturbation theory is applied to a 
Schwarzschild black hole have indicated that an impor- 
tant contribution to the energy radiated dur ing an infall 
has its origin just outside the event horizon [l8|. It is on 
this region we will focus. Since we are interested in the 
plunging phase of the orbit, we make use of the near hori- 
zon approximation, which is introduced in Sec[TTl Fea- 
tures of the final stage of the plunge that influence the 
emitted radiation are identified in Sec lIIII and ApplB] 
The Teukolsky formalism is used to compute the grav- 
itational waveform from the end of the plunging stage 
in Sec lIIIi and to compute the energy and momentum 
fluxes in Secs lIVI and |V] respectively. Numerical explo- 
rations on a Kerr background within the perturbation 
framework have indicated that spin enhances the radia- 
tion emission [l^,[2^. This was understood in the context 
of the lower dampin g ra tes of the quasi-normal modes in 
the presence of spin[l3|. In our work, the spin enhance- 
ment effect is carefully quantifled analytically during the 
end of the plunge phase. 

Our work aims at analytically quantifying features of 
the transition from the plunge to the ring-down phase 
which cannot be captured by extending purely post- 



Newtonian results [2l| . The analysis presented here does 
not extend into the regime of free quasi-normal ringing 
observed when a highly distorted black hole approaches 
Kerr in the absence of an external source. This is due 
to the fact that the perturbing second black hole is still 
present and therefore a source of waves. 

A possible concern about modeling the end of the 
black hole plunge in the extreme mass ratio limit with 
first order perturbation theory is that the problem has 
been oversimplified by ignoring non-linear perturbations 
and employing the point particle approximation for the 
smaller black hole. However, it should be observed that 
the dynamical timescale associated with the small black 
hole is comparatively short, making the equilibration pro- 
cess of internal perturbations of the small black hole rapid 
and thus justifying the point particle approximation. A 
huge advantage of the model is that it takes into account 
the fully relativistic frame dragging effect on the orbit in 
the vicinity of the event horizon. A further advantage 
is that it allows for the exploration of the maximally 
spinning black hole limit. This is a case of particular 
theoretical interest which can only be explored by taking 
a limit within this analytic perturbation framework. Fi- 
nally, the perturbation formalism implicitly takes into ac- 
count the horizon deformation of the background space- 
time, although care has to be exercised when dealing with 
a source that crosses the event horizon, as is shown in 
AppO It is hoped that with the rapid advance of nu- 
merical simulations into the extreme mass ratio regime, 
the predictions of our model will soon be numerically 
tested. 

The analysis presented in this paper allows us to iden- 
tify the following effects. At the end of the plunging 
stage just before the small black hole passes through the 
event horizon, the radiation is dominated by a peculiar 
feature indicative of the ergoregion. This feature is best 
described as follows. In its own reference frame, the small 
black hole passes through the event horizon in finite time; 
however, to a distant observer it remains a radiating mi- 
rage stuck on the event horizon and dragged into mo- 
tion at a tremendous speed, emitting gravitational waves. 
In practice, it is not the mirage emitting gravitational 
waves, but the nonlinear perturbed gravitational field 
that has stored up energy and angular momentum and is 
slowly leaking them out to the observer at infinity. The 
frequency and damping rate display universal behavior 
in that the earlier details of the orbital trajectory are ir- 
relevant; only the position at which the small black hole 
crosses the horizon matters. The emitted radiation has 
a frequency and damping rate analogous to quasi-normal 
ringing but at a frequency distinct from it. Pretorius [3| 
gave a heuristic description of the sinusoidal bobbing mo- 
tion observed in the numerical simulations of two black 
holes in the so-called super kick configurations in terms 
of a frame dragging phenomenon. He mentions that the 
position in the orbit where the merger takes place ulti- 
mately determines the magnitude and direction of the 
kick. While the magnitude of the kick cannot be accu- 
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rately accounted for by the extreme mass ratio analysis 
of the plunge phase of the orbit, the dependence of the 
radiated linear momentum on the angular position on the 
horizon at which the small black hole enters can, and is 
calculated in Sec|Vl We conclude the paper with possible 
observational and analytic implications in Scc lVIl 

Throughout this paper, geometrized units, such as 
G — c = 1 are used. 



II. ORBITAL DESCRIPTION OF PLUNGING 
OBJECT 

In this section we describe the nature of the orbit in the 
final plunging phase. The rapid timescalc on which the 
dynamics of the plunging phase occurs, implies that the 
orbital evolution of the small black hole is well approx- 
imated by a geodesic in the background Kerr geometry, 
ignoring the radiation reaction effect. Since our focus is 
on the late plunging phase, we then further specialize the 
orbit to the near horizon region, also referred to as the 
near-horizon limit. 

The geodesic equations describing the trajectory of a 
freely falling test particle in the Kerr geometry can be 
expressed in Boyer-Lindquist coordinates {t, r, 9, 4>} using 
the equations 
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In the above expressions, the constants of motion for a 
particular orbit are denoted E, and C for energy, 
azimuthal angular momentum and Carter constant re- 
spectively; furthermore, the function A is defined to be 
A = r'^ — 2AIr+a^ and the afhne parameter s is related to 
the proper time r by s = / dr/S, with = r'^+a^ cos^ 6. 

To facilitate understanding of the radiation content 
emanating from the orbit just before the test mass falls 
into the horizon, we further specialize these equations to 
the near horizon limit. Denote the radial position of the 
event horizon hy r = = M + \J M"^ — a? . The quan- 
titative features of an in-falling orbit in the near horizon 



region (r « r+), can be adequately described by approx- 
imating the geodesic equations (|2.ip - (|2.2p by 
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where VLh — 2Mr+ horizon's angular velocity, 

and 6*0 is the polar angle at which the particle falls 
into the horizon, and the dimcnsionless constant k = 
v/1 - {a/M f. ^ 

The geodesic in the near-horizon limit has the following 
analytic solution; 
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where E^nd 0o are integration constants. This approx- 
imation holds for t — to > > r+ . 

Note that the well-known gravitational time dilation 
effect is clearly manifest in terms of the singular behavior 
of Eq. (|2.7p . While the particle passes through the event 
horizon in finite proper time in its own local reference 
frame, it appears to a distant observer who describes the 
particle motion in terms of coordinate time to be slowing 
down infinitely as it approaches the event horizon. As 
a result, the motion in the radial and polar directions 
becomes frozen onto the event horizon at their entering 
positions. If the central black hole has no spin, this parti- 
cle will cease to be a source of gravitational waves during 
the final stages of the plunge to the order of approxima- 
tions made. 

If the central black hole is spinning, however, Ea. (|2.10[) 
indicates that the motion in the azimuthal angle does 
not slow down. Instead the frame dragging effect "pulls" 
the particle into motion around the black hole with an 
azimuthal angular velocity that approaches the horizon 
angular velocity. The particle motion around a spinning 
black hole in the final stage of the plunge is thus, from the 
perspective of a distant observer, entirely dominated by 
this azimuthal motion which is insensitive to the previous 
details of the orbital evolution and characterized only by 
the mass and spin of the central object. This effect is 
highly relativistic, it cannot be captured by an analysis 
based on a post-Newtonian expansion, and is the cause 
of strong gravitational wave emission during the plunging 
phase of a spinning black hole. This statement will be 
made explicit in Sec lIIII 
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FIG. 1: ISCO radius, orbital energy (E) and azimuthal angular momentum (L^) as a function of the spin coefficient q = a/M 



While the constants such as Lz and C describing or- 
bital motion do not enter directly into the description of 
the final stages of the plunging orbit itself, they do make 
their appearance in the dominant terms of the perturba- 
tion theory (see App|^ and are thus required for explicit 
evaluation. To this end, we have selected three plung- 
ing orbits for explicit evaluation: the co- and counter- 
rotating plunging trajectories from the innermost stable 
circular orbit (ISCO) in the equatorial plane (i.e. C = 
and ^0 = ■"'/S), and the orbit plunging from the polar 
ISCO with vanishing z-component of the angular mo- 
mentum, (i.e. Lz — 0)- The constants of motion charac- 
terizing these orbits can be calculated by writing down 
the condition that the orbit remain circular, namely by 
requiring that R{r,E,Lz,C) = dR{r, E, Lz,C)/dr = 
at the ISCO radius r = rjsco- The additional con- 
dition that the orbit remain stable adds the equation 
cPR{r, E, Lz,C)/dr^ = 0. If one further assumes that 
the orbit is cither polar or equatorial, namely that C = 
or Lz — 0, then the three equations allow the three un- 
knowns, E, Lz and riscOi to be expressed in terms of 
the black hole mass M and spin a. 

Fig [H depicts the ISCO radius rjsco, orbital energy 
E and azimuthal angular momentum Lz as functions of 
the spin coefficient q — a/M. Due to spin-orbit cou- 
pling, the orbital energy of the ISCO depends strongly 
on the spin coefficient when the orbit lies in the equato- 
rial plane. When the orbit is co-rotating/counter-rotating, 
the spin-orbit coupling works as a repulsive/attractive 
force that stabilizes/destabilizes the orbit with respect 
to the gravitational attraction of the background black 
hole. As a result, the radius of the co-rotating/counter- 
rotating ISCO becomes smaller/larger when compared 
to the ISCO in the non-spinning limit. This effect is dis- 
played in the first panel of FiglTJ For the polar orbit 
which corresponds to the super kick configuration, the 
orbital energy of ISCO depends weakly on the spin co- 
efficient. It turns out that only the E and Lz orbital 
parameters ultimately enter the waveform calculation in 
Sec lIIIl However for completeness, the Carter constant 



C in the case of the plunge from the polar ISCO orbit 
can be found in terms of energy E and rjsco to be 



III. GRAVITATIONAL WAVES EMITTED 
DURING THE FINAL STAGE OF THE PLUNGE 

The gravitational waves emitted at infinity can be 
computed using the black hole perturbation theory ex- 
pressed using the Teukolsky formalism. The radiation 
contained in the resulting waves is expressed in terms of 
the Newman-Penrose curvature scalar ^'4: 

*4(r->oo) = i (/i+-i/ix) , (3.1) 

where x represents the wave polarization and the 
double dot indicates two time-derivatives. 

Initially it may appear that, due to its bewildering 
complexity, the full waveform calculation of a plunging 
orbit within the context of black hole perturbation the- 
ory may best be attempted numerically by means of di- 
rect integration of the Teukolsky equation as was done 
in [ll, [l^, [23, ■ However, it turns out that identifying 
the dominant terms and spin effects in the calculation an- 
alytically is particularly insightful. Furthermore, in the 
rapidly spinning case, it is highly likely that numerics will 
not be able to capture the final stages of the plunge and 
that this regime is thus only accessible via analytic per- 
turbation theory. The full Teukolsky machinery is sum- 
marized in gory detail in App|^ This section highlights 
the aspects that make the problem tractable analytically, 
the approximations made, and the results obtained. 

The scalar ^"4 can be expanded in terms of its Fourier- 
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Harmonic components as 



tion (13.311 can be evaluated to be 



1 



(r — iacos^ 



(3.2) 



where Rernui{r) and Sirnu>{()) are the radial and angular 
Teukolsky functions respectively. The influence of the 
motion of the point particle as it plunges into the black 
hole enters the radial Teukolsky equation as a source 
term. The particular solution to the radial Teukolsky 
equation can be found by means of evaluating an integral 
containing a Green's function constructed from an expan- 
sion of homogeneous solutions to the Teukolsky equation. 
The fact that we are only interested in the solution at 
r — > cxD simplifies the expression given in ApplAlto 
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where B^"^^ is a mode dependent constant, T^mw the 
source term related to the stress energy tensor describ- 
ing the particle motion and R}^^ is the homogeneous 
in-going Teukolsky function with asymptotic behavior 
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In the above equation, k = uj — ■ma/2Mr^ and r* is the 
tortoise coordinate defined by {dr* /dr) — (r^ + a^)/A. 

The nature of the problem, with the point particle on 
the final stages of the plunging orbit, requires that this 
source integral need only be evaluated near the horizon to 
capture the essentials features of the the radiation emit- 
ted during this phase. Thus the asymptotic form of R}^-^^ 
given in Ea. (|3.4p can be used. It is useful to observe that 
in this region A w 2KM{r — r_|_). 

The dominant contribution to the source term results 
the rapid frame dragging motion observed when the par- 
ticle is frozen onto the horizon at the angle of impact 9q. 
This contribution is contained in the TjTjm component of 
the stress energy tensor projected along the tetrad legs, 
defined in App[Xl The dominant scalings of the other 
components of the source term on the event horizon for 
a plunging orbit are written out in full in AppIbI 

With these approximations, the function Zfmuj enter- 
ing the asymptotic form of the the radial Teukolsky func- 
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where e// = (r — r_|_)/r_|_. 

Using the solution to the geodesic equations special- 
ized to the region near the event horizon (|2.8p - (|2.10p to 
evaluate along the orbit and r* — > [r^/ K)h\(r — r-(-), 
the t-integral of (|3.5p becomes 
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where the lower bound of the integration, T, must be 
chosen so that 6"^'^^''+^'"^"*''' << 1, thus, ensuring the 
validity of the near-horizon expansion. 

From Eq. (|3.7p it becomes clear that the emitted radia- 
tion is strongly peaked at the frequency lu = mQ.H — i-:^- 
Taking this into account, the asymptotic form of the cur- 
vature scalar ^4 is 



'^4^{r — > 00) = — TT 



t*0 



(3.8) 



where the sum is evaluated at an angular frequency of 



H 



'2r+ 



This result represents a decaying 
spectrum of gravitational waves with frequency around 
LD = mVLu being damped in a manner that admits a qual- 
ity factor of 



Q 



(3.9) 



The analytic results characterizing the gravitational 
waves from the plunging phase of the orbit, yields the 
following observations. Firstly, the frequency and qual- 
ity factors of the observed radiation differ from the well- 
known quasi-normal ringing modes of black hole pertur- 
bation theory. The quasi-normal ringing that is observed 
and expected to follow most binary black hole collisions 
can be characterized as damping modes associated with 
vacuum perturbations. The quasi-normal ringing is of- 
ten associated with the relaxation of perturbations of 
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the stationary black hole horizon. The origin of the ra- 
diation described by Eq. (|3.8p is somewhat different. It 
is produced by perturbations coherently induced by the 
particle still present on the horizon. The damping be- 
havior is primarily due to the fact that the particle falls 
through the horizon and ceases to be a source of per- 
turbation. For a quantitative comparison with the least 
damped {I — m — 2) quasi-normal mode, one can make 
used of a fitting formula suggested by Echeverria [3] that 
gives the quality factor and frequency of this mode as a 
function of spin to be 



i^QN = [1 - 0.63(1 



a 



(3.10) 
(3.11) 



where /g and fa are functions of order unity as the spin 
factor ranges from to 1, i.e. fq e [1.05,0.95] and 
fa e [1.02,0.97]. Note that for spin factor a/M < 0.97, 
both the frequency and quality factor of the quasi-normal 
mode exceed those predicted by the driven motion of the 
plunging phase. As a result, the moment the particle 
ceases to be observable and thus to serve as a source 
driving the radiation, the quasi-normal ringing effect will 
rapidly begin to dominate the observables. In the region 
of very high spin (a/M > 0.97) the two effects of quasi- 
normal ringing and radiation driven by particle motion 
on the horizon become almost indistinguishable. In the 
low spin limit the difference between the frequency and 
quality factor of the quasi-normal ringing and those as- 
sociated with the radiation during plunge phase begin to 
differ considerably. However, it should also be noted that 
in this limit, our assumption that the dominant term in 
the stress energy tensor expansion is the frame dragging 
induced rotation breaks down. 

Secondly, while the small particle is plunging, the ra- 
diation emitted broadcasts features of the background 
spacetime, namely Kerr, rather than features reminiscent 
of its orbital trajectory. The spectrum peak at around 
Lu ^ ttiVLh is determined entirely by the constants de- 
scribing the Kerr black hole without any reference to the 
orbital constants of the plunge. As argued in SecHIl the 
frame-dragging effect accelerates the rotational velocity 
of the particle around the spin axis of the black hole. 
The origin of the universal value of the frequency peak 
can thus best be explained as the frame dragging effect 
compensating for the gravitational time dilation effect of 
the radiation. 

Thirdly, the quality factor associated with the radia- 
tion emitted during the plunging phase approaches in- 
finity in the maximally spinning limit. The radiation 
emitted during the plunge in this limit may provide a 
sensitive means of estimating the spin of a rapidly ro- 
tating black hole. The quality factors for the radiation 
of gravitational waves with (m — 2) for a series of spin 
parameters are tabulated below. 



The large quality factor obtained for the rapidly spinning 
limit can be used to explain the enhanced emission of 
gravitational energy and momentum from the plunge ob- 
served in the numerical studies made within the Teukol- 
sky perturbation framework [l^, [23] . The feature that 
the quality factor approaches infinity in the maximally 
spinning limit may appear to endanger the assumptions 
made earlier regarding the insignificance of radiation re- 
action on the orbit. For if a mode at a finite frequency 
lasts an infinitely long time, it will radiate an infinite 
amount of energy. It should be noted that the quasi- 
normal modes of Kerr display a similar feature. In the 
past, it was feared that this infinite quality factor in the 
case of quasi-normal modes implied the instability of the 
Kerr black hole, however this turns out not to be the case 
p^ . [T§ ] . The resolution of this conundrum with respect 
to the current problem is explored in the Sec lIVI 

Fourthly, the relevant expansion parameter near 
the horizon of the black hole is = {''" ~ 

r+)/r+. In the analysis performed in this section, 
only the leading order terms were retained. This 
approach is now justified by the following argument. 
Suppose one goes one step further, and computes 
the higher order terms of the orbit and stress en- 
ergy. This will give rise to additional terms pro- 
portional to J°°dt e*"*-''""«*e-"('=/''+)(*-*")e*'=(*-*") cx 
u-mnH+^n{n/2r+) ^^ich enter the calculation for the ra- 
dial Teukolsky function, Eq. (|3.5p . where n = 2;3,---. 
This suggests that the curvature perturbation in the 
plunging phase can more generally be written as 



£m,n=l,2. 



(3.12) 



a/M 


0.95 


0.97 


0.98 


0.99 


0.995 


0.998 


0.9999 


Q 


6.08 


7.98 


9.85 


14.04 


19.92 


31.58 


141.41 



where the frequency lj entering the expression is different 
for every mode and given by a; = mQn — *'^2^' ■'^'^te 
that in the above equation, the damping factor has been 
enhanced by a factor of n and so the dynamics described 
thus far are truly the dominant behavior, since all higher 
order terms are more strongly damped. 



IV. ENERGY FLUX 

The picture presented thus far is that of a small 
body falling into the central black hole on a dynamical 
timescale with respect to proper time, which for suffi- 
ciently extreme mass ratios, is assumed to be well short 
of the radiation reaction timescale. This justifies the as- 
sumption that gravitational radiation backaction on the 
orbit is negligible and that the particle plunges on a 
geodesic trajectory. 

In Eq. (|3.9p of Scc|TTl it was shown that for rapidly spin- 
ning black holes, the quality factor becomes very large in 
the maximally spinning limit. This fact could possibly 
indicate that the gravitational wave emission at the hori- 
zon angular frequency for a maximally rotating black hole 



7 



may be considerable and last for a long period of time. 
If this were the case, it would nullify the assumptions 
previously made regarding the geodesic trajectory of the 
plunging black hole. 

In this section, we calculate the gravitational energy 
flux during the final stage of the plunge and estimate the 
maximum amount of energy the particle can emit during 
the plunging phase. In doing so, we check the validity 
of the assumptions made and explore the possibility of 
direct detection of this stage of the orbit. 

The gravitational energy flux at infinity can be ex- 
pressed as 



dE 









2' 


lim 


— / dn 


J dm4 




r — >OQ 


16ttJ 







(4.1) 



This suggests that even though gravitational waves could 
be emitted for an extremely long time at a finite fre- 
quency in the maximally spinning limit, the wave ampli- 
tude of this radiation is suppressed and the resultant to- 
tal energy radiated small. This shows conclusively that if 
the mass ratio is sufficiently extreme, the backaction due 
to the emission of gravitational waves is negligible. The 
approximation of the spacetime geometry by the linear 
metric perturbation theory remains valid, as does the as- 
sumption that the in-falling particle moves on a geodesic 
orbit. This result further implies that it will be difficult 
to observe highly relativistic dynamics of the particle or- 
biting just outside the horizon in the maximally spinning 
case because of the strong suppression of the wave am- 
plitude. 

The observed suppression of the wave amplitude in the 
maximally spinning case can be explained using the ra- 
dial Teukolsky equation (|A5[) . Near the horizon the ho- 
mogeneous radial Teukolsky equation (jASP becomes 



where ajm = mflH — i-^^ is used. The total energy radi- 
ated from t — r* — T to t — r* — oo is given by 

I ~ ^ 

S K ^ " 



-e 



16 r+ 



aE- 



E 



|'S'£ma„(6'o)r 



{E-nnL^r 

2 



■ sin-* do 



Tjtrans 



TDin 



(4.2) 



Eq. ()4.2|) can be evaluated by using the analytic expres- 
sions for the homogeneous Teukolsky functions in the low 
frequency limit given in AppO It can be observed that 
the relative contribution of higher £-modes converge very 
rapidly; in fact, they fall off super-exponentially, allowing 
us to consider only the £ = 2 modes. So doing we have 



AE 



1024^ rl 



-e 



{E-QhL, 



sm t/Q 



X < — + sin^(6lo) , , 

I 45 V 15 5 ^ V \kMJ 



\kMJ 

^-8sin^((?o)+ 5sin4(f?o; 



\km) 



/512 56 . ^ 4 . 4,^ s\ / a \M 
+ [^-Y--'iOo) + -s.n^ieo)J (— ) |. 

(4.3) 

In the maximally spinning limit, the scaling behavior of 
the total radiated energy with respect to the quality fac- 
tor can be expressed as 



AE cx cx 



(4.4) 



d / 1 
(r- r+)— - 
dr \r 



dRimui\ , [kr+Y + 2iK{kr+) 



r+ dr 



Rlmuj — . 

(4.5) 



The pole in the potential has a coefficient of . As 
a result, in the maximally rotating limit (k — > 0) it be- 
comes increasingly difficult for radiation originating from 
a source near the horizon to propagate outward and reach 
infinity. The potential barrier in the Teukolsky equation 
thus produces a suppression effect that competes with 
the high quality factor of emitted radiation described in 
Secin] and ultimately dominates as the black hole be- 
comes maximally spinning. The competition between 
these two effects result in a spin coefficient less than unity 
at which the maximum amount of energy is radiated dur- 
ing a plunge. This feature is present in all orbits explored. 

Fig[2] displays the total energy radiated during the 
plunging stage of orbits originating from ISCO. The fac- 
tor e^(''/''+)^-^~*") present in Eq. (|4.3p is weakly depen- 
dent on the orbital trajectory and set to unity for sim- 
plicity. Fig[2] plot (a) shows the energy radiated if the 
original orbit is equatorial. The non-zero z-component 
of the angular momentum (L^) results in finite energy 
being radiated in the limit of zero spin, i.e. a/M = 0. 
As the background spin increases from zero, the radiated 
energy of the co-rotating plunging orbit decreases. This 
effect is due to an effective decrease in the angular mo- 
mentum resulting from the frame-dragging of the back- 
ground black hole - the factor of [aE — L^Y in Ea. (|4.3p 
makes this apparent. In the counter-rotating case, the 
same effect results in a net increase in energy radiated as 
spin increases from zero. As the spin further increases, 
the competition between the increasing quality factor as- 
sociated with the radiation and the suppression effect due 
to the effective potential barrier, described in the previ- 
ous paragraph, sets in. A maximum in radiated energy is 
achieved at a/M = 0.7 and a/M = 0.72 for the co- and 
counter-rotating cases respectively. 
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FIG. 2: Total energy radiated during plunging phase of orbit from ISCO as a function of spin parameter. Plot (a) displays 
the energy radiated when plunging from an equatorial ISCO (C=0) in the co- and counter-rotating cases. Plot (b) shows the 
radiated energy if the particle originates from a polar orbit for various horizon impact angles, ^ = f i f > § 



FiglUplot (b) shows the energy radiated from the end 
of the plunging phase when the particle plunges from a 
polar ISCO orbit. The azimuthal angular momentum 
(Lz) vanishes for this case and the resultant energy radi- 
ated approaches zero in the limit of zero spin, a/M = 0. 
As we increase the spin of the background black hole, 
an effective angular momentum represented by the fac- 
tor {aE — Lz)'^ in Ea. (|4.3p comes into play and the ra- 
diated energy increases. For particles plunging from a 
polar orbit, the radiation suppression effect noted in the 
equatorial case is also observed for high spin factors and 
a maximum in radiated energy is reached. This peak oc- 
curs at a spin factor of a/M = 0.87 when the impact 
angle with the horizon is Oq = or Oq = 7r/3 and at 
a/M = 0.86 when the impact angle is = t^/'^- The best 
fit to the slopes of the three energy profiles in FigOplot 
(b) for the range of spin factor between 0.7 < a/M < 0.8 
is found to be AE = 3.22 x 10"<'(a/M) - 18.94 x 10"^ 
for 6*0 = 7r/6, AE = 1.24 x 10-^{a/M) ~ 7.20 x 10"^ for 
6*0 = 7r/3 and AE = 0.87 x 10-^{a/M) - 5.89 x 10"^ for 
9„ = 7r/3 . 

In the case of a particle plunging from a polar orbit, 
the radiated energy is dependent on the angle made with 
the spin axis when the particle impacts the horizon. The 
angular dependence of the radiated energy is displayed 
in Fig [3] for various values of spin parameter. 

Observe that the maximum in radiated energy does not 
occur when the particle falls in on the equatorial plane, 
^0 = 7i'/2, but that it occurs instead for impact angles 
of 9o — 1.07 and Oq — 2.06. One would naively expect 
to find that the radiated energy is largest in the config- 
uration where the particle has maximum orbital velocity 
V when it impacts the horizon, namely on the equatorial 
plane at = 7r/2. An initial scaling argument would be 
that the energy radiated is roughly proportional to the 
square of the metric perturbation which in turn is pro- 
portional to the square of the stress energy tensor of the 




0.5 1 1.5 2 2.5 3 



Plunging angle 9 

FIG. 3: Total energy radiated during plunging phase of a 
particle originating on a polar ISCO as a function of impact 
angle with the horizon. Radiation for background Kerr black 
holes with spin parameter (a/M) = 0.8, 0.9, 0.99 are plotted. 

source term T ozv^ (x (wr)^ sin^ 6*0. This effect is indeed 
displayed by the appearance of the factor of sin^ 6q in 
Ea. (j4.3p . The energy emission is, however, complicated 
by the angular pattern resulting from the spin-2 nature 
of gravitational waves, which is described by spheroidal 
harmonics that enter Ea. (|4.3p . The origin of the double 
maxima in the emitted energy can be explained by focus- 
ing on the (^TO)-sum of Eq. (|4.2p . Making use of equation 
([Ug|) . to evaluate the ratio B*™"V B™, the sum can be 
expressed as 

Y,\0Jrnf\Slmoj'', (4.6) 

tm 

with Re[u)rn\ = mfln- The dominant contribution to the 
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sum in Eq. (|4.2p and thus the radiated energy is made 
by modes with the largest |m| value for a given I. If 
we consider the contribution from (I — 2)-modes using 
the low-frequency result in App|Cl the 0o-dependence be- 
comes 

sin4(0o) {(1 - cos(0o))' + (1 + cos(0o))'} , 

which has two maxima at 6*0 = 1.06 and Bq = 2.08 and a 
local minima at = 7r/2, and thus accurately explains 
the ^o-dependence observed in Fig[31 



V. LINEAR MOMENTUM FLUX 

During the evolution of a typical binary system, the 
longest period of time is spent in the inspiral phase. It 
is during this phase that most of the gravitational en- 
ergy is emitted. By comparison to the inspiral phase 
the energy flux, computed in Sec IIVI emitted during the 
plunging phase is insignificant. The situation is very dif- 
ferent when one considers the linear momentum flux. The 
net linear momentum emission during the inspiral phase 
is small, because the binary orbit gradually progresses 
through a series of quasi-stable circular orbits resulting 
in nearly isotropic momentum emission. As a result, the 
momentum emitted during the plunging stage is impor- 



tant and may be comparable to that emitted during in- 
spiralling phase [ig] . In this section we quantify the lin- 
ear momentum flux carried away by gravitational waves 
during the final plunge phase of a spinning black hole. 
The emitted linear momentum flux is given by 



dt 



lim 

r — >oo 



— - / dn n., 




2" 


1 dt^i 




167r J 


J —oo 





(5.1) 



where nj = (sin 6 cos 4>, sin 9 sin 0, cos 9) . 

The z-component of the momentum flux, dPz/dt, van- 
ishes near the horizon. This is yet another manifestation 
of the phenomenon discussed in SecITTl where the particle 
is frozen onto the horizon due to the time dilation effect 
and its motion is dominated by the velocity component 
in the ^-direction. In the near horizon expansion of the 
source term given in ApplB] at the leading order only, 
the particle motion in the x- and y-directions are taken 
into account and the net emission of z-linear momentum 
vanishes as a result. 



The calculated total linear momentum carried away by gravitational waves during the period T < t ~ r* < cx3 is 

/ d9.in\9)S,„..^{9){S,„,.^,{9)r , 

where we set m' — to + 1. To facilitate explicit evaluation, the analytic expressions for the Teukolsky functions in the 
low frequency given in AppICjare used. The resulting expression is 

MPx + ^Pv) ~ M V r:^ sm^ 9q 

f 16 .16 / a \ /16 8 . 2rn^\ ( V 

" 1 135 - ^45 l^dVf ) + U - 15 l^J 
. / 176 16 . ^^ / a \3 



-i-i-^(^o))(^)V(li-isin,.))(^ 



(5.2) 



The momentum radiated during the flnal plunging energy of the same orbit, Fig[2] (a), the linear momentum 
phase of an orbit originating at ISCO on the equatorial plotted here displays a finite value in the limit of a non- 
plane is plotted in Fig |4] (a) . Analogous to the radiated spinning Kerr background hole due to the orbital angu- 
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FIG. 4: Total momentum radiated during the final plunging phase of orbit at ISCO as a function of spin parameter. Plot 
(a) displays the momentum radiated when plunging from an equatorial ISCO (C — 0) in the co- and counter-rotating cases. 
Plot (b) shows the radiated momentum if the particle originates from a polar orbit (Lz = 0) for various horizon impact angles, 

n 7r TT 7r 

6 ' 3 ' 2 



lar momentum of the plunging particle. The momentum 
radiation further displays features reminiscent of the ra- 
diated energy, in that the radiated momentum decreases 
for the co-rotating plunge and displays an increase at first 
followed by a decrease for the counter-rotating plunge as 
the background spin increases from zero. A feature dis- 
tinct from those displayed by the radiated energy in Fig. 
[2] (a) is that no local maxima of the radiated momentum 
is apparent. This is due to the factor 1/(1 — i{a/2KM)) 
which comes from the time integration of the momentum 
flux. This can be better understood by observing that 
in Ea. (|5.ip . the momentum flux displays damped oscilla- 
tory behavior and for large spin, the oscillation frequency 
becomes high. When taking the integral to obtain the to- 
tal radiated momentum, the positive and negative parts 
of the integrand cancel, resulting in a small value for the 
integral. This stands in contrast to the radiated energy 
where the integrand is a positive decaying function with- 
out any oscillatory behavior. 

The oscillatory behavior of the momentum flux made 
explicit in Eq. (|5.ip indicates that the linear momentum 
flux is emitted in various directions in x/y-plane. This 
feature is to be expected if the particle is in effect pulled 
into orbit around the black hole as it plunges and so 
the radiated momentum has an oscillatory nature that 
will reflect the horizon velocity of the flnal black hole. A 
similar process was alluded to in Fig. 14 in the second ref- 
erence of to explain the anti-kick phenomena observed 
in [13. 

The momentum radiated during the flnal stages of 
plunge from a polar ISCO is shown in In Figd] (b). 
Analogous to the radiated energy plotted in Figl2] (b), 
the radiated momentum increases from zero as the back- 
ground spin increases. The radiated momentum reaches 
a maximum at a/M = 0.82, when the in-falling angle is 



^0 = ■"'/S or 6q = tt/3 and at a/Af — 0.81 when the in- 
falling angle is Oq = tt/2. As the spin-factor ranges over 
0.55 < a/M < 0.75, the radiated momentum can be flt- 
ted by AP = 1.23x 10-^(a/M)-5.65x 10"** for 6*0 = n/6 
AP = 0.95 X 10-^(a/M) - 4.12 x lO""^ for 9o = n/3 and 
AP ^ 0.94 X 10-^{a/M) - 3.62 x lO"'^ for 6*0 = 7r/3. 




FIG. 5: Total momentum radiated during plunging phase of 
a particle originating on a polar ISCO as a function of impact 
angle with the horizon. Radiation for background Kerr black 
holes with spin parameter [a/M) — 0.8, 0.9, 0.99 are plotted. 

The dependence of the radiated momentum on the an- 
gular position at which the particle impacts the horizon 
during a plunge from a polar orbit is show in In FiglS] 
The angular profiles of the radiated momentum (Figl5|) 
and those of the radiated energy (Figl3|) differ slightly. 
The features of the local minimum at Oq — and the 
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two maxima off the equatorial plane are robust. For the 
case of radiated momenta, the maxima occur at = 1-33 
and ^0 = 1-82. Following an argument similar to that em- 
ployed in the case of the radiated energy, the observed be- 
havior can be explained by analyzing the 0o-dependence 
of the {££'m)-sum in the first line of equation (15. 2|) . The 
dominant functional behavior is found to be 

sin5(0o) {(1 - cos(0o))' + (1 + cos(0o))'} , 

which admits maxima at 0o — 1.35 and 6o = 1.79 and a 
local minima at — 7r/2. 

VI. CONCLUSION 

This paper explores the dynamics of the final stages of 
the plunge of a small black hole into a Kerr black hole 
and characterizes the resulting gravitational radiation, 
using perturbation theory. The purpose of this explo- 
ration is to gain analytic understanding of the features 
which dominate the dynamics of this transition phase of 
the merger process before the onset of quasi-normal ring- 
ing and to identify any universal features. Our analysis 
strives to highlight the effects of the spin of the central 
black hole on the resulting gravitational waveform and 
momentum flux while building an intuitive picture of the 
processes responsible for their origin. It should be noted 
that the all graphs plotted in this paper and comments 
made are relevant only to radiation emanating from the 
so-called near horizon region during the plunging phase. 

In Sec lIIIl it was found that the waveform resulting 
from the plunging phase can be expressed as 

h{t)^ ^ V.,„)e-*"<'"-'* , (6.1) 

n,m— 1,2--- 

where the frequencies entering the expansion take on the 
universal form uj(,n,n) — fnilHt — inK/(2r^). Therefore, 
the time dependence of the emitted radiation is entirely 
determined by the characteristics of the massive central 
black hole. This waveform describes the transition pe- 
riod between radiation dominated by the inspiral and 
radiation dominated by the free quasi-normal ringing. 
The inspiral stage broadcasts the characteristic frequen- 
cies of the orbiting particle and is strongly affected by 
the constants of motion describing the orbit. The quasi- 
normal ringing phase is characterized by the frequencies 
and damping rates associated with free vacuum pertur- 
bations of an excited Kerr black hole as it settles down to 
equilibrium. The frequencies that enter expression (j6.ip 
can be ascribed to neither of the two preceding scenar- 
ios. They describe a state of forced oscillation distinct 
in frequency from free quasi-normal ringing but are at 
the same time devoid of any characteristic identifying 
the source of the perturbation. In fact, the source can 
be best described as a reluctant mirage on the horizon 
- a mirage dragged into motion and then retained for 
distant observation by an all-dominating frame dragging 



effect. The frequency associated with originates 
from the azimuthal motion of the orbit on the horizon, 
which is entirely determined by parameters of the back- 
ground black hole. The method by which this mirage is 
stripped of its identity while it in essence passes through 
the event horizon and is assimilated into the large black 
hole is described in mathematical detail in SecHIl 

The wave amplitude /i(m,n) still retains characteristics 
of the source and depends on the orbital constants. The 
damping rate associated with W(m.n) indicates that the 
radiation corresponds to a decaying part of the wave- 
form. A feature that this frequency shares with the quasi- 
normal mode frequencies is the existence of a high quality 
factor in the maximally spinning limit a/M 1, for ex- 
ample, Q = 14.04 for a/M = 0.99 and Q = 141.41 for 
a/M = 0.99. A more detailed comparison between the 
least damped quasi-normal mode frequency and u>(^2,i) is 
given toward the end of Sec lIIIl 

The gravitational flux of energy and linear momen- 
tum originating from the plunging phase is computed 
in Sec lIVI and Sec|Vl respectively. Orbits plunging from 
ISCO in the equatorial plain (i.e. C = 0) and from a 
polar ISCO (i.e. = 0) have been considered. In gen- 
eral, the following characteristics can be identified. The 
rapid increase in quality factor is counteracted by a term 
in the amplitude /i(m,ri) that rapidly approaches zero as 
spin approaches maximality. This ensures that all forms 
of radiation approach zero in the limit of maximal spin 
and that the total emitted energy and momentum dur- 
ing this phase remain finite. This effect further ensures 
that a maximum in emitted energy as a function of spin 
is attained in all cases considered. For the plunge from 
an equatorial orbit the maximum in radiated energy oc- 
curs around a/M — 0.7 and for the plunge from polar 
orbit the maximum occurs for a spin parameter around 
a/M = 0.86. 

The plunge from a polar ISCO orbit is of particular 
interest because it can be used as a toy problem to illu- 
minate some of the features observed in numerical simu- 
lations of black hole mergers in super-kick configurations. 
In particular, the dependence of the emitted linear mo- 
mentum on the angle 9q made with the rotation axis when 
the particle enters the horizon is computed in this paper. 
Contrary to expectation, it was found in Sec|V] that the 
produced profile does not have a simple sinusoidal depen- 
dence on angle, but rather that it is complicated by the 
spin-2 nature of gravitational radiation, as displayed in 
Figli 

The radiated linear momentum originating from a po- 
lar plunge reaches a maximum with respect to the spin 
parameter at around a/M = 0.86. While the origin of 
this peak is well understood within the context of our 
model, the question remains as to, whether the peak is 
an artifact of linear perturbation theory and our underly- 
ing assumptions or whether one will find its counterpart 
in numerical simulations. We also find that the radi- 
ated energy and linear momentum can be fitted by linear 
functions of a/M in the domain 0.7 < a /A! < 0.8 and 
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0.55 < a/M < 0.75, respectively. 

A final observation regarding the computed momen- 
tum flux of SeclVlEq. llS.ip is the existence of oscillatory 
behavior with a frequency we now consider characteris- 
tic of a plunging event involving spin. The introduced 
phase implies that the linear momentum flux is emitted 
in various directions in the x/y-plane, consistent with our 
picture of a miraged particle stuck on the event horizon 
radiating momentum in a manner that reflects the ve- 
locity of the final black hole while it slowly recedes from 
view. This could provide a heuristic explanation of the 
origin of the anti-kick phenomenon observed in |17| and 
is suggestive that the anti-kick is indicative of the final 
merger event. 
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where the potential V is defined to be 
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and the function K — (r^ -I- a?)u) — ma. The source 
term T^mw contains the effect of the stress-energy tensor 
perturbing the background spacetime. The relationship 
between the source term and the stress-energy tensor is 
given by 
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APPENDIX A: TEUKOLSKY FORMALISM 

This appendix summarizes the Teukolsky formalism 
used in the derivation of the Weyl scalar \E'4 which con- 
tains the radiation content of the perturbed spacetime 
with a prescribed source term. 

The scalar ^'4 can be decomposed into its Fourier- 
harmonic components as follows 



*4 = 



1 



(r. 



lacosf 
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da;e'''^'+'^'^Ri„Ur)SimUO)- 
(Al) 



where Simuj (0) and Remui (r) represent the angular and ra- 
dial Teukolsky functions respectively. The angular func- 
tion, obeys the angular Teukolsky equation for s = —2, 
explicitly. 
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where the potential U is given by 

(m -2cosef 
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U = a^uj^ sin^ I 



sm a 
2 — 2maLU - 
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(A3) 



and A denotes the eigenvalue of Sgmu. The angular 
Teukolsky function is normalized so that. 



sin 6 de\Si. 
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(A4) 



Ttm^ = 4 J dndtp-^p-\Ti + T2)e- 
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(A7) 



2V2 



(A8) 



T, = --p^pA'V+ip-^V+ip-^pTrr 



(A9) 

where we use p = {r — iacos9)~^ ^ p — (r + iacos9)~^ , 
and the operators I?+ = dr + iK/A , and Cg — dg + -^-g — 
aLusinO + scot 9. The terms r„„, Tmn and Tmm denote 
the tetrad components of the stress-energy tensor. In 
this paper, a point particle was used as the source of the 
metric perturbation. The stress-energy tensor associated 
with a point particle can be expressed as 



dtY^ dzf" dz" 
S sin 9 Kdr J dr dr 



X S{r - r{t))5{9 - 9{t))5{cf) ~ cl){t)) , (AlO) 



where p is the mass of the particle. Using the geodesic 
equations (|2.ip - (|2.2p . the source term can be simplified 
to, 



Time. = M / dte'"*-™^A2 [dl {A,,^2<5(r- - r{t))} 

+ dr {{Afanl + ^mml)'5(r - r{t))} 

)(5(r - r(<))] 



(All) 
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where 
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(A16) 

(A17) 

(A18) 

(A19) 
(A20) 



Here the operator £j is defined by £| ^ dg — + 
acu sm9+s cot 9 , indices have been dropped from the an- 
gular Teukolsky function, i.e. S = Simuj ^md t ~ dt/dr = 
{l/Y){dt/ds). 



For a given source the particular solution of the radial 
Teukolsky equation can be found by means of a Green's 
function method. To this end, define two homogeneous 
solutions to the radial Teukolsky equation which satisfy 
the boundary conditions 
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for r ^ cx) ' 

(A22) 



in these expressions k — uj — maj^Mr^ and r* is the 
tortoise coordinate defined by {dr* /dr) = (r^ + a2)/A. 
In terms of homogeneous solutions the particular solution 
of the radial Teukolsky function can be expressed as 
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—RZ^ / drRZ^A-'T,^^ , (A23) 



where the Wronskian Wimuj is given by 

Wimu = 2iw-B^;;^C|™^'' . Of particular interest for 

the purposes of gravitational wave extraction is the 
nature of the solution at infinity, 
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In the event that the source term is a point particle mov- 
ing along a geodesic, such as the scenario considered in 
this paper, the function Zi^uj can be expressed as 
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(A25) 

(A26) 

(A27) 
(A28) 



the terms (|A27I) and (|A28[) appear as the boundary terms when equation (jA24[) is integrated by parts. If the source 
particle orbit is bound within the radial domain, namely somewhere between infinity and the horizon, these terms 
vanish. 

For the source under consideration in the present paper, care has to be exercised since the orbit crosses the horizon 
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and the term (jA28|) may contribute to the integral. It is however found that this term too vanishes, which is reasonable 
since no waves propagate from the horizon out to infinity. 



APPENDIX B: LEADING ORDER BEHAVIOR OF SOURCE TERMS 



The source terms generated by a particle plunging on a geodesic orbit near the horizon of a black hole display the 
following leading order behavior when expanded with respect to distance from the horizon 



^nnO 
^mm2 



0((r-r+)3) 
o^kEisco 



SVC? 9q 



4Mr4 



0(r 



(r+ — ia cos 

^mnO 



+ 0((r-r+n 



KO^Eisco I kr 

0{r^r+), 

. KO^Eisco kr+ r+ 




sin^ OqS- 



r — 



+ 0{{r-T+f)., 



2\f^Mr\ K r++ia cos do 
Ko^Eisco — io. cos 



la Costco • 7 n ^ ^/ ^ 
sm-^ 6*05^ + 0{r - r+) , 



4V27rMri r+ + ia cosOq 



sm 



•oS{r-r+)+0{{r-r+f), 



(Bl) 
(B2) 
(B3) 

(B4) 
(B5) 
(B6) 

(B7) 



Note that the dominant contribution of the particle's geodesic motion enters the source terms by means of 
This term dominates due to the rapid azimuthal velocity in duced by the frame dragging effect. 



r 



APPENDIX C: TEUKOLSKY FUNCTIONS IN 
THE LOW-FREQUENCY LIMIT 

This appendix summarizes the the analytic considera- 
tions required to compute the functional form of homo- 
geneous Teukolsky functions in the low-frequency limit. 



1. Angular Teukolsky Functions 

The solutions to the angular Teukolsky equation are 
also known as spin-weighted angular spheroidal func- 
tions. The analytic structure of these functions was stud- 
ied in [2^. By expanding the spin- weighted spheroidal 
functions in terms of Jacobi polynomials the angular 
Teukolsky equation can be reduced to three-term re- 
currence relations that are given by Eqs. (21-24) of [23j . 
The solutions of the recurrence relations were shown, in 
Eq.(25) of 23], to converge very rapidly with respect to 
an expansion in frequency, i.e. with respect to expansion 
parameter auj. As a result, in the low-frequency limit, 
Jacobi polynomials accurately capture the analytic be- 
havior of the spin- weighted spheroidal harmonics. 

To leading order in the low frequency limit, the angu- 
lar Teukolsky functions can thus be well approximated by 
spin-weighted spherical harmonics. It is this approxima- 
tion that is adopted in this paper to aid computational 
simplicity. In particular we have, 



where f > |m| and i = 2,3, ■ ■ ■ . The low £, spin- weighted 
spherical harmonics used in the computation are now 
given 



^2Y2±2iO,0)^y^il±COSe)\ 



(C2) 



-2i^2±i(6',0) = iw-sin6»(l±cos6'), (C3) 

4 V TT 
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-2>2o(0,O) = -Y— sm- 



(C4) 



2Yim{0,O), ii-m + 2), (CI) 



2. Radial Teukolsky Functions 

The properties of the solutions to the radial Teukolsky 
equation used to aid calculation are briefly summarized 
in this section. The analysis of the analytic structure 
of the radial Teukolsky equation was first performe d by 
Leaver [23| and further developed by Mano et. al .[25j. 
it is this reference to which the reader should refer for a 
treatment more detailed than the synopsis given here. 

The approach adopted in the analysis is to expand the 
homogeneous solutions in terms of two different series ex- 
pansions, each strongly convergent at a particular bound- 
ary condition. Near the horizon hypergeometric func- 
tions were used (refer to Sec. 2 of Ref. [25|) to capture the 
behavior of in-going and outgoing waves and implement 
the boundary condition on the horizon. An expansion in 
Coulomb functions was then used to give an asymptotic 
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expression at infinity, and to fix the boundary condition 
of no in-going waves (refer to Sec. 3 of Ref. |25|). These 
two expansions where then matched in the interior to 
obtain a solution vahd over the whole domain obeying 
the appropriate boundary conditions on both boundaries 
(refer to Sec.4 of Ref. [111). 

Following an approach similar to that employed in the 
case of the angular Teukolsky function, the radial Teukol- 
sky equation can be reduced to three-term recurrence 
relations valid for both expansions. Once again the so- 
lution to the three term recurrence relations shows very 
rapid convergence with respect to a low frequency expan- 
sion. In particular the expansion parameter e = 2ujAI is 
introduced. 

In the calculation performed in this paper, the ratio of 
the incoming wave amplitude, B"^ to the amplitude of 
the transmitted wave, iJ*'''^"* as defined in Eq. (jA2ip is 
required. Using [2^ the ratio of these wave amplitudes 
can be obtained in terms of a power series expansion in 
e. The rapid convergence of the expansions employed in 
[25I I with respect to e allows us to consider only the lead- 
ing order term in the expansion to make a quantitative 
estimate of the required ratio. The wave amplitude ratio 



so calculated is 

j^trans 



(i; + 2)!(^-2)!r(^+l + »mf) 
(2£)!(2£+l)! r(3 + zm^) ^ ^ ' 

(C5) 



where q = a/M and k = -^/l ~ q^. The factorials in the 
denominator of Eq. (jC5[) . cause the computed ratio to 
decrease very rapidly as i increases. As a result only the 
lower ^ contributions need be considered. Specifically, 
the lowest two terms contributing to the expression are 
computed to be 



^trans 



-I— w5e-*^-'"'"(l + 0(e)) 



B' 



1=2 



i=3 



15 
2 

'945^ 



(C6) 



'(Ke)e-^^ « (3 + tm^^ (1 + 0(e)) 

(C7) 



[1] Frans Pretorius, Relativistic Objects in Compact Bina- 
ries: From Birth to Coalescence, Colpi et al.. Springer 
Verlag, Canopus Publishing Limited, jarXiv:0710. 1338 1 
and references therein 

[2] M. Campanelli, C. O. Lousto, Y. Zlochower, D. Mer- 
ritt Astrophys.J. 659 (2007) L5, J. A. Gonzalez, M. D. 
Hanna,, U. Sperhake, B. Brugmann and S. Husa, Phys. 
Rev. Lett. 98, 231101 (2007), M. Campanelli, C. O. 
Lousto, Y. Zlochower, D. Merritt Phys. Rev. Lett 98, 
231102 (2007). 

[3] J. G. Baker, W. D. Boggs, J. Centrella, B. J. Kelly, S. T. 

McWilliams, M. C. Miller, J. R. van Meter Astrophys. J. 

682, L29 (2008) 
[4] S. Dain, C. O. Lousto, Y. Zlochower Phys. Rev. D 78, 

024039 (2008) 

[5] J. D. Bekenstein J D Astrophys. J. Lett. 183, 657 (1973). 
[6] M. J. Fitchett, Mon. Not. R. Astron. Soc. 203, 1049 

(1983), A. G. Wiseman Phys. Rev. D 46, 1517 - 1539 

(1992). 

[7] L. E. Kidder, Phys. Rev. D 52, 821 (1995), S. H. Miller, 
R.A. Matzner arXiv:0807.3028, to be publish at Gen. Rel. 
Grav. 

[8] L. Blanchet, M. S. S. Qusailah, C. M. Will Astrophys.J. 

635 (2005) 508, T. Damour, A. Gopakumar Phys.Rev. 

D73 (2006) 124006. 
[9] J. D. Schnittman, A. Buonanno Astrophys. J. Lett. 662 

L63 (2007), J. D. Schnittman, A. Bounanno, J. R. van 

Meter, J. G. Baker, W. D. Boggs, J. Centrella, B. J. 

Kelly and S. T. McWiUiams submitted to Phys. Rev. D, 

[arXiv :0707.030f 
[10] M. J. Fitchett and S. D. Detweiler Mon. Not. R. Astron. 

Soc. 211, 933 (1984), M. Favata, S. A. Hughes and D. E. 

Holz Astrophys. J. Lett. 607, L5 (2004). 
[11] L. Blanchet, Living Rev. Relativity 9, (2006), 4. 

|http:/ /www.livingreviews.org/lrr-2006-4^ and references 

therein 



[12] M. Sasaki and H. Tagoshi, Living Rev. Relativity 6, 
(2003), 6. http://www.livingreviews.org/lrr-2003-6, Y. 
Mino, M. Sasaki, M. Shibata, H. Tagoshi, T. Tanaka 
Prog.Theor.Phys.Suppl. 128 (1997) 1-121, and references 
therein 

[13] S. Chandrasekhar and S. L. Detweiler, Proc. R. Soc. Lon- 
don A344, 441(1975), S. L. Detweiler, Proc. R. Soc. Lon- 
don A352, 381(1977), S. L. Detweiler, Astrophys. J. 239, 
292 (1980), E. W. Leaver, Proc. R. Soc. London A402, 
285(1986). 

[14] F. Echeverria, Phys. Rev. D40, 3194 (1989). 
[15] V. Ferrari and B. Mashhoon, Phys. Rev. Lett. 52, 1361 
(1984) 

[16] H. P. NoUert, Phys. Rev. D47, 5253 (1993) 
[17] M. Koppitz, D. PoUney, C. Reisswig, L. RezzoUa, J. 
Thornburg, P. Diener, E. Schnetter Phys. Rev. Lett. 99, 
041102 (2007), D. PoUney, C. Reisswig, L. RezzoUa, B. 
Szilagyi, M. Ansorg, B. Deris, P. Diener, E. N. Dorband, 
M. Koppitz, A. Nagar, E. Schnetter Phys. Rev. D 76, 
124002 (2007) 

[18] M. Davies, R. Ruffini, W. H. Press and R. H. Press, Phys. 

Rev. Lett. 27, 1466 (1971), M. Davies, R. Ruffini and J. 

Tiomno, Phys. Rev. D5, 2932 (1972). 
[19] T. Nakamura and M. Sasaki, Phys. Lett. A 89, 185 

(1982) , T. Nakamura and M. Haugan, Astrophys. J. 269, 
292 (1982), T. Nakamura and M. Sasaki, Prog. Thore. 
Phys. 67, 1778 (1982). 

[20] Y. Kojima and T. Nakamura Phys. Lett. A 96, 335 

(1983) , Y. Kojima and T. Nakamura Prog. Theore. Phys. 
71, 79 (1984). 

[21] M. Sasaki Prog. Theor. Phys. 69, 815 (1983). 
[22] K. Oohara and T. Nakamura Prog. Theore. Phys. 70, 757 
(1983). 

[23] E. D. FackereU and R. G. Grossman, J. Math. Phys. 18, 
1849 (1977). 

[24] E. W. Leaver, J. Math. Phys. 27, 1238 (1986). 



16 



[25] S. Mano, M. Suzuki and E. Takasugi, Prog. Theor. Phys. 
95, 1079 (1996). 



